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ABSTRACT 


The  simple  H&ckel  method  was  first  applied  to  the  electronic 
structure  of  the  iron-porphyrin  complexes  by  Pullman  et  al.  0  In  this 
paper,  their  work  is  extended  to  include  (a)  the  effect  of  a  dipole  or  a 
point  charge  placed  at  the  sixth  coordination  position,  and  (b)  the  effect 
of  a  nitrogen  atom  placed  at  the  fifth  coordination  position.  A  set  of  new 
parameter  values  is  used,  whose  estimation  is  made  by  directing  special 
attention  to  their  dependence  on  the  charge  distribution  among  the  atoms. 

The  resulting  charge  distribution  for  ferro-porphyrin  seems  to 
be  reasonable.  The  fact  that  the  position  of  the  Soret  peak  is  insensitive 
to  the  sixth  ligand  can  be  understood  from  the  resulting  orbital  energy 
levels. 

The  difficulty  of  finding  a  reasonable  charge  distribution  for  ferri- 
-porphyrin  is  discussed. 


1.  INTRODUCTION 

One  of  the  molt  interesting  chemical  properties  of  re -porphyrin 
complexes  ie  the  existence  of  "high-spin”  and  "low-spin"  states.  For 
example,  the  number  of  unpaired  spins  in  Hb  -  HjO  and  Hi  -  H20  is  4  and 
5  respectively,  while  that  in  Hb  -  02  and  Hi  -  OH*  is  0  and  1  respectively. 

Another  interesting  feature  is  their  spectra.  All  these  complexes 
show  an  extremely  strong  hand,  called  the  Soret  band,  at  390  *  430  mp.  The 
position  of  the  peak  does  not  change  very  much  for  Hb  -  I^O,  Hb  -  02  and 
even  for  metal -free  porphins. 

The  complexes  are  known  to  play  very  important  roles  in  oxidative 
and  reductive  reactions  in  biological  systems.  It  has  been  pointed  out  that 
this  redox  power  is  closely  related  to  the  orbital  energies  of  the  highest 
occupied  and  the  lowest  unoccupied  molecular  orbitals.  ^ 

It  is  very  difficult  to  give  a  theoretical  interpretation  of  these  phenom¬ 
ena  from  first  principles.  And  it  is  even  more  difficult  to  explain  the  oxygen- 
-uptake  power  of  haemoglobin.  However,  it  is  clear  that  one  of  the  moot 
fundamental  pieces  of  knowledge  about  this  molecule  is  its  electronic  struc¬ 
ture. 

We  use  the  simple  Hilckel  method  in  this  work  not  because  it  is  ade¬ 
quate  for  our  purposes,  but  because  it  is  simple  enough  to  be  easily  applied 
to  complicated  systems  such  as  the  iron -porphyrin  complexes. 

The  simple  Hdckel  method  was  first  applied  to  the  iron-porphyrin 
complexes  by  Pullman,  Spanjaard  and  Berthier.  ^  Their  calculations  are 
extended  in  this  work  to  include  (a)  the  effect  of  a  dipole  or  a  point  charge 
placed  at  the  sixth  coordination  position,  and  (b)  the  effect  of  a  nitrogen  atom 
placed  at  the  fifth  coordination  position.  We  also  introduce  new  parameter 
values,  whose  estimation  is  made  by  directing  special  attention  to  their 
dependence  on  the  charge  distribution  among  the  atoms. 

In  Pullman  and  others'  work,  ^  it  was  assumed  that  in  the  high- 
-spin  complexes  the  3d  electrons  of  the  iron  atom  do  not  participate  in  the 
bonding  between  Fe  and  the  porphin  ring,  while  in  the  low  -spin  complexes, 
they  do.  The  calculated  molecular  orbital  energies  differed  considerably  for 
these  two  types  of  complexes.  The  work,  therefore,  tug gested  that  when  a 
complex  changes  its  spin  state,  there  are  drastic  changes  both  in  the  charac- 
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ter  and  in  the  anargy  of  the  molacular  orbital*  involving  the  3d  orbital*. 

Thi*  may  wall  ba  trus,  bat  the  specific  aaaumption  mentioned  above  is  rather 
artificial.  Hare  we  lift  this  assumption.  Changes  in  molecular  orbitals  will 
be  brought  about  through  the  dependence  of  the  parameters  on  the  charge 
distributions  which  are  different  for  the  high*  and  low-spin  state. 


2.  THE  MODEL 

All  the  peripheral  substituents  of  the  porphyrin  ring  are  neglected  so 
that  the  system  under  consideration  contains  20  C's,  4  N's,  and  one  Fe  in 
the  molecular  plane  (cf.  Fig.  1).  and  one  N  under  the  Fe  atom  -  this  wiil 
be  denoted  hereafter  N'  •  The  assumed  distances  between  the  atoms  are 
the  following:  Fe  -  N  1.80  A,  N  -  C'  1.38  A,  C'  -  C  and  C'  -  C" 

1.4  A  ,  Fe  -  N'  2. 10  A .  The  bond  angle  ZC'NC'  of  108°  is  also  assumed. 
These  figures  are  fairly  close  to  the  observed  ones  for  nickel  ethioporphyrin 
and  myoglobin.  ^ 

The  ligand  at  the  sixth  coordination  position  apparently  determines 
the  spin  state  of  the  system.  We  replace  this  ligand  by  a  dipole  or  a  point 
charge  and  investigate  its  effect.  A  dipole  with  its  axis  perpendicular  to  the 
molecular  plane  or  a  point  charge  is  thus  placed  above  the  Fe  atom  at  the 
distance  of  2.0b  A  .  The  system  has  the  symmetry  C^y  . 

The  system  is  so  large  that  only  the  bonds  between  Fe  and  four 

surrounding  N's  are  considered  for  the  v  part;  also  the  possible  v  bonding 

between  Fe  and  N'  is  neglected.  With  all  foese  simplifications,  there  are 

42  electrons  (one  each  from  C,  three  each  from  N,  8  from  Fe  and  2  from  N') 

for  the  ferrous  complex  and  41  electrons  for  the  ferric  complex.  To  describe 

them,  we  use  the  following  38  atomic  orbitals:  one  2pv  of  each  C,  one  2p* 

and  one  sp^  hybridised  <r  of  each  N,  five  3d,  one  4s  and  three  4p  of  Fe  and  one 
2 

sp  hybridised  orbital  of  N'  pointing  to  Fe  . 
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3.  METHOD  OF  CALCULATION 

Our  model  has  the  symmetry  of  C^y  but  if  we  disregard  the  point 
charge  or  dipole  at  the  sixth  coordination  position  and  the  N*  atom  at  the 
fifth  coordination  position,  the  symmetry  is  °4h-  It  is  more  convenient  to 
name  some  of  the  molecular  orbitals  after  the  irreducible  representations  of 
the  group  The  names  of  the  irreducible  representations  of  the  groups 

and  C^y  and  the  relation  between  them  can  be  seen  in  Table  I. 

TABLE  I  Irreducible  representations  of  the  groups  and  C4y  • 


Putting  aside  the  orbital  of  N'  ,  we  form  linear  combinations  of  the 
37  atomic  orbitals  available'  The  linear  combinations  form  bases  for  the 
irreducible  representations  of  .  These  are  called  symmetry  orbitals 

and  explicit  forms  are  given  in  Appendix  I.  The  irreducible  representations 
which  appear  and  number  of  linearly  independent  functions  belonging  to  these 
follow:  Alu(2),  A2u(5).  Blu(3),  BZu  (3)  ,  E  (7)  .  A  (3).  B  (2).  B2  (1) 
and  Eu  (2)  .  *  * 

When  the  atom  N’  and/or  a  point  charge  or  a  point  dipole  are  con- 
sidered,  we  have  Aj(9)  instead  of  Aj  (3)  and  A2u  (5)  and  an  sp 
hybrid  of  N'  pointing  to  the  Fe  atom.  Off-diagonal  elements  between 
Bjg  and  B2u  ,  B2g  and  B  Ju  ,  and  Eu  will  be  aero  under  the  approx¬ 
imations  employed  in  the  present  calculations.  (This  will  be  explained  in 
the  next  section).  Even  in  the  general  case  we  shall  keep  the  symmetry 
notation  of  except  for  A  and  A2u  . 

Thus  we  have  expressed  the  symmetry  orbitals  in  terms  of  atomic 
o&itals.  In  order  to  know  how  much  these  symmetry  orbitals  mix  to  give 
molecular  orbitals,  we  have  to  set  up  and  solve  the  secular  equations. 

Integrals  which  appear  in  the  secular  equations  can  be  reduced  to 
those  involving  atomic  orbitals.  In  the  simple  H&ckel  theory,  they  are  the 
Coulomb  integrals  a  ,  the  resonance  integrals  p  and  the  overlap  integrals 
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S  •  They  are  defined  &• 


CC  p 

f  X  p  (  u 

Heft  (  u  Xp  i  i ; 

(3. 1) 

Pn  * 

J  X  y  c  n 

Hett  1  1  X-  i  1  J 

jVl  , 

(3.2) 

S  p<|  -* 

W“' 

X i  i )  Jvt# 

(3.3) 

where  Xp  and  Xq  are  the  atomic  orbitale.  The  total  Hamiltonian  ia  assumed 
to  be  a  sum  of  one -electron  effective  Hamiltonians  H£££  ; 

H  -  Z.  t  i )  .  (3.4) 

i 


The  one -electron  effective  Hamiltonian  (i)  may  be  considered 
to  have  the  following  form: 


H  .ft  c  i  )  -  -  A i  h  v*  c  i  )  ' 


(3.5) 


where  is  th.  effective  potential  due  to  atom  a  . 

The  Coulomb  integral  ap  for  an  atomic  orbital  Xp  centred  on  the 

atom  a  can  be  split  into  two  parts,  a  '  and  a  "  : 

P  P 

OCp  »»  j  X  p  (  *  j  A  v  Lfj.  )  X  p  dv  +2.  j  Xp  vTj,  X  p  <A  >f 

^  (3.6) 

“  0Cp  *  OCp,"  . 


The  first  term  a  '  may  be  equated  to  a  suitable  valence-state  ionisation 

‘  lj 

potential  for  occupied  atomic  orbitals  (e.g.  2pw  of  C  ,  N  and  3d  of  Fe  ) 
or  to  an  electron  affinity  for  empty  atomic  orbitals  (e.g  4s  and  4p  of  Fe++)  . 
The  second  term  ap”  is  neglected  unlesis  atom  b  has  an  appreciable  formal 
charge.  When  the  atom  b  does  have  a  formal  charge,  Zb  ,  the  term  is 
evaluated  by  assuming  that  a  point  charge,  Zb  ,  is  situated  at  the  nucleus  b  . 

The  resonance  integrals  8,  are  estimated  by  a  slightly  modified 

Ph 

Wolfsberg-Helmhols  approximation,  •  namely 


Pm 
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t  K 
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(3.7) 


•5- 


where  K  ■  2  for  ir  bond*  and  1. 67  for  v  bond*.  The  baste  of  this  approx¬ 
imation  is  not  quite  clear  but  it  is  certainly  most  convenient  in  reducing  the 
number  of  unknown  parameters.  The  modification  we  made  consists  in 
separating  off  the  electrostatic  part  from  a^'.  This  is  made  in  order 
to  fulfill  the  requirement  that  the  resulting  orbital-energy  differences 
should  not  change,  if  we  change  the  sero  point  of  the  energy  by  adding  a  con¬ 
stant  potential  to  the  Hamiltonian. 

The  values  of  the  overlap  integrals  can  be  evaluated  without  too 
much  difficulty,  at  least  in  principle. 

It  should  be  noted  that  all  the  parameters  a  ,  p  and  S  depend  on 
the  formal  charge  which  the  atoms  carry.  This  necessitates  the  use  of  the 
method  of  trial  and  error  which  goes  as  follows.  We  first  assume  a  charge 
distribution  and  determine  the  parameters.  Then  we  obtain  the  molecular 
orbitals  and  orbital  energies,  and  thus  we  are  led  to  a  plausible  electron  con¬ 
figuration.  Analyzing  the  wave  function  of  this  electron  configuration,  we 
obtain  a  resulting  charge  distribution.  If  this  is  close  to  the  starting  charge 
distribution,  we  can  stop  there.  If  the  resulting  electron  distribution  differs 
from  the  starting  one,  we  have  to  start  again  assuming  a  more  appropriate 
charge  distribution.  This  is  a  requirement  of  self-consistency  in  charge 
distribution  which  seems  to  be  quite  important  in  dealing  with  very  hetero- 
-polar  systems  like  metal -porphyrine  complexes. 

The  potentials  due  to  a  point  charge  (chargeZe),  a  dipole  (moment  p 
with  the  axis  on  the  z  axis)  and  an  axial  quadrupole  (moment  Q  *  with  the 
axis  perpendicular  to  the  z  axis  and  parallel  to  the  y  axis)  placed  at  the 
distance  R  from  the  origin  can  be  expressed,  respectively,  in  the  polar 
coordinates  as  follows; 
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The  quadrupole  moment  is  defined  as 
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where  p  is  a  charge  distribution  function. 
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In  the  last  expression  (9  ,  <p  )  is  defined  as 


In  calculating  matrix  elements  of  these  potentials  with  respect  to  the 
s-  ,  p-  and  d-type  atomic  orbitals,  the  integrals  can  be  carried  out  over  the 
angular  parts  easily  and  the  resulting  expressions  are  listed  in  Appendix  II. 


4.  VALUES  OF  THE  PARAMETERS  FOR  FERROUS  COMPLEXES 

As  was  explained  in  the  previous  section,  our  parameter  values  depend 
on  the  charge  distribution  we  assume. 

The  choice  of  parameters  by  Spanjaard  and  Berthier  ^  indicated  that 
they  started  from  a  model  in  which  all  the  atoms  are  neutral.  The  resulting 
formal-charge  distributions  were  Fe^’®+  porphin**®’  and  Fe1,5+  porphin1*® 
for  low-  and  high-spin  ferrous  complexes.  In  view  of  these  results,  we 
start  with  a  model  in  which  the  formal-charge  distribution  is  FeZ*®+  N®’ 5'  , 
all  carbon  atoms  being  considered  as  neutral. 

4.1.  Valence -state  ionisation  potential  and  electron  affinity 

The  first  terra  ap’  of  the  parameter  u  (cf.  (3.6))  is  equated  to  a 
suitable  valence -state  ionisation  potential  for  occupied  atomic  orbitals  or  to 
an  electron  affinity  for  unoccupied  atomic  orbitals.  The  values  of  these  used 
in  the  present  calculation  are  collected  in  Table  11. 


TABLE  II 


Valence -state  ionisation  potentials  and  electron  affinities 
(in  eV). 


Ionisation  potential  I 

Electron  affinity  A 

Atom 

Orbital 

I 

Atom 

Orbital 

A 

l  .. 

C 

2pv 

11.42  ^ 

Fe2+ 

4s 

16. 18  3) 

N' 

2pr 

14.26  2* 

Fe2+ 

4P 

11.41  ^ 

N' 

s 

18.67  2) 

n°.5- 

2ps 

7.79  2* 

, 

n°.5- 

r 

12.08  2) 

Fe2+ 

3d 

30.65 

3 

The  valence  state  is  assumed  to  be  the  tetra-valent  sp  state.  The 
value  is  taken  from  the  paper  by  H.O.  Pritchard  and  H.A.  Skinner, 
Chem.  Rev.  55,  745  (1955). 

For  nitrogen,  the  valence  state  is  assumed  to  be  a  mixture  of 

s2p3  with  factor  i/3  and  V,  sp4  with  factor  it/3*  1°  estimating  I  for 

2  J  2  3 

the  sp  hybridised  r ,  I  of  2s  for  V,  s  p  is  assumed  to  be  the  same 

4  3 

as  that  for  Vj  sp  .  Numerical  values  are  taken  from  Pritchard  and 
Skinner's  paper. 

For  N®’  an  additional  assumption  is  necessary.  It  is 
observed  that  the  ionisation  potential  of  2p  electrons  for  N"  ,  N  and 
N+  is  a  linear  function  of  the  ionicity  of  the  atom.  Thus  it  is  assumed 
that 

I(N0-5’).£tI(N)  +  I(N")}  * 

I(Fe2+  3d)  »  E(3d5)  -  E(3d*)  , 

A(Fe2+  4s)  •  £(34*)  -  £(3d64s)  , 

A(Fe2+  4p)  ■  E(3d*)  -  E(3d*4p)  . 

Atomic  term  values  are  found  in  the  book  by  C.E.  Moore;  "Atomic 
Energy  Levels",  National  Bureau  of  Standards,  Vol.  II  (1952). 
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4.2.  Contribution  to  a  from  the  other  etome 
carrying  a  formal  charge 


The  second  term  ie  estimated  by  assuming  the  point  charges  Ze 
(Z  *  +2  for  Fe,  and  -0.  5  for  each  of  four  N)  are  located  on  the  nuclei.  For 
C  .  for  example,  we  have 


-  2 


(  »  2p*e  (f)!1 

J  1  r  -  Rfcf 

(  » 

J  I  P  -  R'ni 


4  ^ 
dtr 


r  i»»'c  «»r 

J  If-  R„l 


d  tf  — 

(4. 1) 


where  R  m  and  R'n  are  the  position  vectors  of  the  nitrogen  nuclei  nearer 
and  farther  to  the  C  atom  under  consideration,  respectively.  Each  term  in 
the  right  hand  side  of  (4.  l)  is  a  so-called  nuclear  attraction  integral. 


For  the  atomic  orbitals  used  in  estimating  these  integrals,  for 
example  2pv^,  in  (4.  1),  we  take  the  Slater  orbitals.  The  orbital  exponents 
6  are  assumed  as  follows: 


6(N)  *  1.95, 

6(C)  a  1.625, 

6  (Fe  3d)  a  2.08, 

6  (Fe  4s)  a  6  (Fe  4p)  =  1. 01 


(4.2) 


For  N  and  C  ,  the  values  of  the  nuclear  attraction  integrals  are  obtained 
by  interpolation  of  the  table  by  Kotani  et  al.  ^  To  calculate  the  nuclear 
attraction  integrals  which  involve  Fe  atomic  orbitals,  the  formulae  given 
in  Appendix  II  are  used,  where  the  values  of  the  incomplete  T  function  are 
interpolated  from  Pearson's  table.  ^ 
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4.3.  Value*  of  Coulomb  Integrals  a. 

The  values  of  a  '  .  a  ”  as  well  as  a_  (»  o'  4  a" )  are  listed  in 
P  P  P  P  P 

Table  I II. 


TABLE  III  The  parameter  a  for  the  model  Fe2+N^' 5’  (in  eV). 


Atom 

Orbital 

a' 

aM 

a  »  o'  4  a" 

n°-5- 

w 

-7.79 

-7.92 

-15.71 

n°-5- 

<r 

-12.08 

-11.72 

-23.80 

C' 

IT 

-11.42 

0.91 

-10.51 

C" 

tr 

-11.42 

0.44 

-10.98 

c 

IT 

-11.42 

0.28 

-11. 14 

Fe2+ 

3d 

-30.65 

15.80 

-14.85 

Fe2+ 

4s 

-16. 18 

13.  12 

-3.06 

FeZ+ 

4p 

-11.41 

13.  12 

1.71 

N' 

-18.67 

-4.99 

-23.66 

4.4.  Overlap  integrals 

The  values  of  the  overlap  integrals  (except  those  involving  N'  orbital) 
were  calculated  by  Spanjaard  and  Berthier  ^  using  Slater  orbitals.  Those 
involving  N'  orbital,  we  have  calculated  using  the  same  AO's. 

We  adopt  Spanjaard  and  Berthier' s  values  for  the  overlap  integrals 
between  C  and  C  ,  and  C  and  N  .  However,  their  values  for  the  overlap 
integrals  between  Fe  and  N  are  modified  by  multiplying  them  by  the  con* 

stant  factor  0.75.  This  factor  is  the  ratio  a/b  .  a  is  the  overlap  integral 

24  91  2 

between  the  Fe  Hartree  -Fock  3d<r  '  and  the  N  sp  hybridised  Slater 

orbital  evaluated  along  the  internudear  axis,  b  is  the  same  integral  be¬ 
tween  the  Fe  Slater  3dr  orbital  and  the  N  sp2  hybridised  Slater  orbital. 
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Thia  procadura  ia  introduced  to  obtain  approximate  overlap  integrala  between 
Fe2+  and  N  ,  but  ia  cartainly  a  vary  cruda  ona.  However,  it  ia  hoped  that 
the  reaulta  of  the  calculation  are  not  ao  aenaitive  to  the  valuaa  of  the  overlap 
integrala.  The  valuaa  uaad  are  given  in  Table  IV. 


TABLE  IV.  Valuea  of  overlap  integrala 

(C  2pw|C2p«) 

0.2444 

(C  2pir|N  2pw) 

0. 1930 

(Fe  3dw|N  2p«) 

0. 1040 

(Fe4pw|N2pw) 

0. 1297 

(Fe  Sdj,  |N  2p<r) 

-0.1133 

(Fe3dxfc.yll|N2pa) 

0. 1963 

(Fe  4a  |N  2p«r) 

0.2129 

(F  e  4pcr  |  N  2p«r) 

0.2939 

(Fe  3d^  |N') 

0. 1560 

(Fe  4a  IN') 

0.2025 

(Fe4p|N') 

0.2976 

4.5.  Re aonance  integrala 

Reaonance  integrala  between  non -neighbour a  are  neglected.  Thoae 
between  neighboura  are  calculated  from  Coulomb  and  overlap  integrala  by 
using  the  formula  (3.7). 
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4. 6.  Effect*  of  a  point  charge  or  a  dipola 

The  cases  in  which  a  dipole  of  either  0.  S  or  1-0  a.u. ,  or  a  point 
charge  of  1.0  a.u.  is  placed  at  the  sixth  coordination  position  are  investi¬ 
gated. 

The  effects  of  these  static  fields  on  the  Coulomb  integrals  (a")  are 
calculated  by  the  formulae  given  in  Appendix  □.  The  resonance  integrals  are 
changed  accordingly  (cf.  (3.7)  ).  However,  direct  contributions  of  the  static 
fields  to  the  resonance  integrals  are  neglected. 

In  the  non-oxygenated  condition  of  haemoglobin,  the  water  molecule 
is  attached  at  the  sixth  coordination  position,  while  in  the  oxygenated  con¬ 
dition,  the  oxygen  molecule  is  attached  there. 

The  permanent  dipole  of  the  water  molecule  is  known  to  be  0.728  a.u. 

(=  1.85  Debye)  by  experiment.  Induced  dipoles  of  the  water  and  oacygen  mole¬ 
cules  placed  at  a  distance  of  2.01  A  from  Fe  are  estimated  to  be  0. 54  and 
0.55  a.u. ,  respectively.  In  this  estimation,  we  must,  of  course,  take  into 
account  the  effect  from  the  negative  formal  charge  of  0.  5  on  each  nitrogen 
atom. 

Assume  that  we  can  neglect  the  polarisation  of  the  surroundings  al¬ 
though  this  is  very  doubtful.  Then  the  main  part  of  the  electrostatic  interaction 
between  the  haem  and  the  oxygen  (or  water)  molecule  may  be  represented 
by  placing  a  dipole  of  0.  5  (1.0)  a.u.  at  the  sixth  coordination  position. 

5.  RESULTS  AND  DISCUSSION  FOR  FREE 
PORPHIN  AND  FERROUS  PORPHYRIN 

5. 1.  Orbital  energy  and  atomic  population  analysis 

Orbital  energies  are  tabulated  in  Table  V,  and  the  significant  energy 
levels  are  illustrated  in  Figure  3.  Because  a  comparison  of  the  absolute 
energies  of  cases  having  different  external  fields  has  little  meaning,  the 
orbital  energy  differences  are  recorded  with  respect  to  the  lalu  orbital,  as 
shown  in  the  figure. 

For  free  porphin,  we  choose  a  model  in  which  the  two  protons  inside 
the  porphin  ring  are  completely  screened  by  tr  electrons.  The  four  nitrogens 
contribute  6  v-electrons  to  the  system.  26  tr-electrons  occupy  the  lowest 
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PABLE  V  Orbital  energies  (in  eV). 


/ 
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13  molecular  orbitals  as  follows: 

( l»2u)2(  ■*,)4  ( <  ■ ">  lu>2  <2',)4  <2‘2«)2  <2b2u)2  <*•/ 

<l*.a)2<3*2»)2  • 

If  we  apply  Muliiken'  ■  population  analysis  to  the  wave  function 
representing  this  electronic  configuration,  we  obtain  the  following  charge 
distribution: 


n1.60  <.0.95  c0.95  <."1.03  . 

The  charge  distribution  suggests  that  the  nitrogen  atoms  attract 
electrons  from  the  carbon  atoms  but  the  degree  is  very  small,  and  the  whole 
molecule  is  almost  homo-polar. 

For  ferrous  porphyrin,  the  degree  of  mixing  of  the  iron  orbitals  in 
the  molecular  orbitals  is  of  some  interest.  In  Table  VI,  the  atomic  popula¬ 
tions  are  given  for  the  molecular  orbitals  which  have  more  than  10  o/o  Fe  3d 
population. 

The  molecular  orbitals  which  have  the  largest  contribution  of  Fe  3d 
atomic  orbitals  are 


2blg 

•V 

3dxt,yi, 

6al 

***»  » 

4e 

g 

*v 

Mxx  ’  3dyz 

and 

lb, 

2g 

♦V 

3d 

xy 

They  are  indicated  so  in  Figure  3.  There  is  one  exception.  In  the  case  of  the 

point  charge  field,  5e  instead  of  4e  has  the  largest  3d  ,  3d  character. 

g  g  xx  ys 

The  order  of  the  orbital  energies  t  of  these  orbitals,  namely 

£(1b*a;  <  8  (4e3)  <  £C6ai>  <  6  <  2  b  ij  )  , 

is  in  agreement  with  the  order  indicated  by  electron  spin  resonance  experi¬ 
ment  on  ferrite  complexes.  1  ^  However,  the  calculated  difference  between 
these  levels  is  much  too  big  in  comparison  with  the  values  obtained  by 
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TABLE  VI  Atomic  population  for  tome  molecular  orbitals. 


j)  No  external  field 


;\Atom 

■■■■ 

j  MOX^ 

Fe(3d) 

Fe(4a,4p) 

4N<r 

4Nw 

4c 

8C' 

8C" 

N‘ 

■V. 

6a. 

0.821 

0.  179 

__ 

• 

• 

• 

• 

0.848 

0.029 

0.069 

V  ~r’ 

0.001 

0.000 

0. 000 

0.053 

0.429 

- 

- 

ym 

0.001 

0.233 

- 

’  2Cg 

0. 114 

• 

- 

0. 177 

0.334 

0.371 

- 

lb, 

2g 

le 

8 

lal 

lbl8 

_ 

1.000 

- 

- 

- 

- 

- 

- 

• 

0.375 

- 

- 

0.579 

0.001 

0.044 

0.001 

- 

-0.008 

Kga 

0.000 

0.000 

0.000 

0.000 

0.307 

0.  179 

“ 

U 

“ 

• 

' 

• 

(ii)  Dipolar  field  of  0.  5  a.u. 


\Atom 

MO^X 

Fe(3d) 

Fe(4s,  4p) 

4N<r 

4Nir 

4C 

8C' 

8C" 

N* 

2b, 

is 

B 

- 

n 

• 

• 

mm 

B 

6ai 

0.032 

i  i 

0.000 

0.000 

0.000 

I 

4c 

g 

B9 

- 

■ 

0.  190 

0.001 

0. 144 

‘■4 

■a 

- 

- 

1 1 

1C8 

mm 

I 

0.616 

1 

1 

0.003 

1 

!  lai 

BBS 

-0.007 

1 

0.000 

0.000 

I 

i  % 

i _ 

0. 170 

■ 

• 

■ 

■ 

- 

B 

TABLE  VI  (coat.) 
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(iii)  Dipolar  field  of  1.0  a.  u. 


IS 

Fe(3d) 

Fe(4t,4p) 

4Nr 

4Nv 

4C 

8C' 

8C" 

N’ 

2blg 

6a  i 

5e 

g 

lb2g 

le 

g 

n 

0. 162 

m 

m 

• 

_ 

BBS 

0.035 

0.052 

0.000 

0.000 

0.000 

0.000 

0.039 

mm 

- 

- 

0. 121 

0.285 

0.153 

0.266 

- 

0.442 

- 

- 

0. 120 

0.013 

0.213 

0.211 

- 

1.000 

- 

- 

- 

- 

- 

- 

- 

0.272 

- 

- 

0.646 

0.004 

0.074 

0.005 

- 

% 

0.  162 

0.838 

(iv)  Point  charge  of  -1.0  a.u. 


Fe(3d) 

Fe(4«,4p) 

4N«r 

4Nir 

4C 

8C’ 

8C" 

N* 

2b. 

!g 

6a, 

0.844 

0. 156 

M 

M 

» 

• 

0.890 

0.033 

0.049 

i 

0.000 

1 

0.000 

0.028 

*g 

4eg 

lb2g 

1Cg 

lblg 

0.522 

- 

- 

0. 154 

0. 133 

0. 165 

- 

- 

0.136 

BE.- 

0.  334 

1.000 

- 

- 

- 

- 

- 

- 

0. 177 

- 

- 

0.627 

0.016 

0. 141 

0.038 

0.156 

• 

0.844 

• 

• 

• 

• 

analysis  of  observed  g  tensor  of  fsrrihaemoglobin  aside.  l2^  For  example, 
the  calculated  value  of  £  (4eg)  -  £  (lb2g)  is  2.5  eV  while  Griffith's 
estixnats  is  only  about  0.2  eV. 

It  is  apparent  in  Figure  3  that  the  orbital  energies  of  the  four  orbitals 
of  strong  Fe  character  are  more  sensitive  to  the  external  field  strength 
than  the  energy  of  the  other  molecular  orbitals.  This  is  quite  natural  since 
the  Fe  atom  is  the  nearest  to  the  source  of  the  external  field.  However, 
both  the  energy  differences  between  these  levels  and  the  percentage  of  Fe 
AO's  in  the  orbitals  do  not  change  very  much  when  the  dipole  is  increased 
from  0  to  1  a.u.  (about  2.  5  Debye)  or  even  for  the  case  when  a  point  charge 
increases  from  0  to  1  a.u. 


5.2.  Electronic  configuration  and  charge  distribution  of  the 
low -spin  state 

For  the  low-epin  state  of  ferro-porphyrin  (with  S  «  0)  ,  the  electrons 
occupy  the  lowest  21  molecular  orbitals  as  follows: 

(lblg)2  Ua/  (l«u)4  <2a/  (leg)4  (3a/  (lb2/ 

(2eg)4  (lb2g)2  (lblu)2  (4a/  (Zb2/  (3eg)4  (5a/ 

(l«lu)2  <«V4  * 

This  configuration  yields,  by  the  population  analysis,  the  charge 
distributions  which  are  given  in  Table  7.  These,  except  for  the  point  charge 
case,  are  reasonably  close  to  our  starting  assumption  about  the  formal 
charge,  namely  Fe2*®*  N®*  ®  .  In  the  case  of  the  point  charge  of  -1  a.u., 
its  effect  is  so  big  that  electrons  are  repelled  from  the  central  Fe  atom  to 
the  peripheral  carbon  atoms,  and  the  formal  charge  on  the  Fe  atom  is  as 
great  as  44.0  .  As  no  stable  ferrous  complex  with  a  negative  ion  at  the 
sixth  coordination  position  is  known,  an  effort  to  make  the  calculation  con¬ 
sistent  in  the  charge  distribution  is  not  attempted. 
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TABLE  VII  Charge  distribution  of  ferro -porphyrin  in  the  low -spin  state. 

The  values  in  parentheses  are  the  formal  charges. 


■ 

No  external  field 

D  =  0.  5  a.u. 

D  * 

1.0  a.u. 

P  «  - 

1.0  a.u. 

B 

1.03 

(-0.03) 

1.05 

(-0.05) 

1.07 

(-0.07) 

1.20 

(-0.20) 

9 

0.88 

(  0.  12) 

0.89 

(  o.u) 

0.92 

(  0.08) 

0.94 

(  0.06) 

H 

1.01 

(-0.01) 

1.02 

(-0.02) 

1.05 

(-0.05) 

1.22 

(-0.22) 

N 

3.66 

t 

O 

• 

O' 

2 

3.64 

(-0.  64) 

3.61 

(-0.61) 

3.53 

(-0.53) 

Fe 

6.16 

(  1-84) 

5.93 

(  2.07) 

5.54 

(  2.46) 

3.97 

(  4.03) 

N’ 

1.91 

(  0.09) 

1.92 

(  0.08) 

1.94 

(  0.06) 

1.97 

(  0.03) 

5.  3.  Electronic  configuration  of  the  high-spin  state 

For  the  high  spin  state  (with  S  =  2),  we  must  have  four  singly -occupied 
orbitals.  The  sum  of  orbital  energies  is  higher  than  that  in  the  low-spin  state, 
and  the  first  stabilizing  factor  we  can  think  of  is  the  exchange  interaction.  Then 
the  singly -occupied  orbitals  should  be  fairly  localized  around  the  Fe  atom.  The 
following  configuration  suggests  itself 

...  (lalu)2(4egy!(6a1)(2blg).  *> 


For  the  case  of  the  point  charge, 

•••  (*«1u)2(Segft^i)(2big) 

would  be  the  configuration,  but  we  shall  not  discuss  it  because  of  its 
inconsistency  in  charge  distribution  with  the  starting  model. 


IS* 


The  charge  distributions  obtained  from  the  above  configuration  are 
given  in  Table  VIII.  The  computed  charge  distribution  ie  not  ae  close  to  the 
starting  charge  distribution  as  it  was  in  the  case  of  the  low-spin  state.  We 
should  repeat  the  calculation  starting  from  a  new  charge  distribution  something 
like  Fe*‘  N®* C°'  *+.  Then  we  could  obtain  new  orbitals  and  their  orbital 
energies  and  see  how  much  they  ere  different  from  those  in  the  low-spin  case. 
We  have  not  repeated  the  calculation,  however,  because  the  dependence  of  the 
molecular  orbitals  and  their  energies  on  the  spin- state  is  caused  also  by 
electron-electron  interaction  which  is  not  explicitly  taken  into  account  in  the 
Hflckel  method. 


TABLE  VDI  Charge  distribution  of  ferro-porphyrin  in  the  high-spin 

state.  The  values  in  parentheses  are  the  formal  charges. 


No  external  field 

D  = 

0. 5  a.u. 

D  = 

1.0  a.u. 

C" 

0.97 

(  0.03) 

1.00 

(  0.00) 

1.02 

(-0.02) 

C' 

0.86 

(  0. 14) 

0.86 

(  0. 14) 

0.86 

(  0. 14) 

C 

1.01 

(-0.01) 

1.02 

(-0.02) 

1.04 

(-0.04) 

N 

3.60 

(-0.60) 

3.60 

(-0.60) 

3.60 

(-0.60) 

Fe 

7.00 

(  1.00) 

6.73 

(  1-27) 

6.40 

(  1.60) 

N' 

1.96 

(  0.04) 

1.97 

(  0.03) 

1.98 

(  0.02) 

The  orbital  energy  differences  (cf.  Table  V)  8  (ba^)  +  6  (2b  ^)  - 
-  2  £(4e^)  are  4.5  eV,  5.2  eV  and  5.9  eV  for  0=0,  0.5,  1.0  a.u. ,  respectively. 
It  seems  not  unreasonable  to  suppose  that  4.5  «v  6  eV  could  be  gained  from  the 
electron  interaction  energy  in  going  from  the  low-spin  to  the  high-spin  state. 

To  settle  this  point,  as  well  as  to  obtain  better  molecular  orbitals,  it  is 
necessary  to  apply  a  method,  such  as  Pariser-Parr-Pople's  in  which 
the  Coulomb  repulsion  between  electrons  is  explicitly  taken  into  account. 

It  may  be  worth  mentioning  the  orbital  energy  difference  £  (6a + 

4  £  (Zbig)  *  2  £  (4%  g)  is  bigger  when  the  dipolar  field  strength  D  is  greater. 

As  there  is  no  big  change  in  the  iron  AO  percentages  of  the  three  molecular 
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orbital*  6a  j,  2b  jg,  and  4*^'  whan  D  Increase*  from  0  to  1  a.u.  (•** 
Table  VI),  there  should  not  be  much  change  in  the  electron  interaction  energy. 

It  looks,  therefore,  like  the  high-spin  state  is  more  unstable  when  the  dipolar 
field  is  greater.  However,  ferrohemoglobin  (with  the  water  molecule  as  the 
sixth  ligand)  is  known  to  be  in  the  high- spin  state  while  oxygenated  hemoglobin 
(with  the  oxygen  molecule,  which  has  only  an  induced  dipole  as  the  ligand) 
is  in  the  tow- spin  state.  Thus  the  computed  result  contradicts  the  experimentally 
known  fact.  This  might  suggest  that  there  is  a  weak  chemical  bond  between 
O2  and  haem  and  that  the  interaction  is  not  essentially  electrostatic. 


5.4.  Spectrum 

Iron-porphyrin  complexes  show  either  one  or  two  bands  of  moderate 
intensity  in  their  visible  spectra  in  the  region  6300  A  (  ~  2.0  eV)  to  5000  A 
(~  2.5  eV)  and  an  extremely  strong  band,  called  the  Soret  band,  at  about 
4000  A  (^  3.0  eV).  Thus  HbOg.  HbH^O  and  even  metal -free  porphyrin  have 
Soret  peaks  around  4000  A. 

The  simple  H&ckel  method  employed  here  is  known  to  be  not  reliable 
in  discussing  spectra.  Nevertheless  we  try  to  see  what  we  can  say  about  the 
spectra  from  the  results  of  the  present  calculations. 

We  assume  that  the  Soret  band  for  metal-free  porphin  is  assigned  to  the 

transitions  from  the  ground  state  to  two  excited  configurations  as  follows:  '  in 

one  an  electron  is  excited  from  3a,  to  4e  ,  in  the  other  from  la.  to 
14)  2u  g  lu 

4eg  ' .  For  iron-porphyrin,  then,  it  is  natural  to  assign  the  Soret  band  to  the 

transitions  5aj -*  5e^  and  laiu-*5e^,  since  5a j,  laJu  and  5e^  correspond 

to  3a2u>  laiu  *eg  in  porphin,  respectively.  We  see  this  not  only  from 

Figure  3,  but  also  from  Table  DC  which  shows  the  atomic  populations  of  these 

orbitals. 

The  differences  between  the  orbital  energies  of  these  two  configurations 
£  (5eg)  -  £  (laiu)  £  (5®g)  "  £(5aj)  ar®  listed  in  Table  X.  This  table 

shows  that  these  orbital  energy  differences  are  insensitive  to  the  dipolar  field 
strength.  There  is  2  suggestion  that  the  inter-configuration  matrix  element  is 


The  effect  of  a  quadrupole  is  examined  and  found  to  be  very  small  for  a 
reasonable  range  for  the  strength  of  the  quadrupole  field. 
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TABLE  IX  Atonic  population*  of  ion*  of  the  highest  occupied  and 

lowest  unoccupied  molecular  orbital*. 


TABLE  X 


(In  eV) 


£  (5«g) ’ t  0*lu)  «<*  2  (5eg)  -  E  (5a  j) 


roughly  constant  for  all  porphins  14\  If  we  assume  this,  Table  X  indicates 
that  the  peak  of  Soret  band  is  insensitive  to  the  external  field  or  the  ligand 
at  the  sixth  coordination  position,  in  accordance  with  the  experimental  observa¬ 
tion. 

The  inter -configuration  mixing  matrix  element  of  0.42  eV  was  used 
by  G  outer  man  14).  This  value  and  the  calculated  orbital  energy  differences 
give  the  peak  of  Soret  band  of  about  2. 5  eV  (~  4960  A),  which  is  too  big  by 
0.5  eV.  However,  the  agreement  is  as  good  as  we  could  expect  in  the  simple 
HO  eke  1  method 


As  the  system  is  large  and  all  the  molecular  orbitals  involved  in  the 
transitions  are  spread  over  the  porphin  ring,  we  need  not  worry  too 
much  about  a  singlet -triplet  splitting.  The  exchange  integrals  between 
these  orbitals  should  be  fairly  small. 
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6.  A  DIFFICULTY  WITH  FERRIC  IRON -PORPHYRIN 

Aa  was  shown  in  Table  VH,  we  obtained  a  charge  distribution  which  is 
reasonably  close  to  the  starting  model  Fe^+  N®‘  ^ '  (or  the  ferrous  iron- 
porphyrin  in  the  low-spin  state* 

We  have  failed  to  obtain  this  consistency  for  the  ferric  iron-porphyrin* 

By  starting  from  the  model  Fe^+  N^'  ,  all  the  orbitals  with  strong  Fe  3d 

character  turn  out  to  have  much  too  low  energiee.  The  resulting  formal  charge 
on  Fe  is  nearly  aero.  By  starting  from  the  model  FeZ+  instead,  the 

orbitals  of  strong  iron  character  still  come  out  much  too  low. 

If  we  regard  the  ferric  complex  as  the  ferrous  complex  minus  an  elec¬ 
tron  in  the  highest  occupied  orbital,  we  arrive  at  the  following  model: 

N  1+  Fe^'  ^  N®*^~  C  *+  (an  electron  missing  from  4e^)  • 

This  model  indicates  that  the  electron  removed  in  forming  the  ferric  complex 
mainly  comes  from  C'  carbons. 

If  we  speculate  about  the  ferric  complexes  in  the  high -spin  state,  the 
above  models  lead  to  a  picture  in  which  unpaired  electrons  are  fairly  delocalised. 
Then  the  problem  immediately  arises}  what  is  the  stabilizing  factor  for  the 
high-spin  state,  as  the  exchange  interaction  should  be  small? 

On  the  other  hand,  if  we  assume  the  unpaired  electrons  are  localised 
around  the  Fe  atom,  we  have  to  explain  why  an  electron  can  be  taken  away 
only  from  the  Fe  atom  (without  disturbing  the  balance  between  the  Fe  atom  and 
the  porphin  ring),  when  the  ferrous  complex  is  oxidized. 

It  would  be  extremely  nice  if  we  could  find  out  the  distribution  of  these 
unpaired  electrons  in  the  complex  from  some  experiment,  such  as  a  double 
resonance  experiment. 
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APPENDIX  I 
Symmetry  Orbitals 

The  co-ordinate  system  and  the  numbering  of  atoms  in  the  molecular 
plane  are  shown  in  Figure  2. 


Orbitals  of  the  v-character 

Alu 

C'  2pw 

1  -  2  +  3  -  4  + 

5  -  6+  7  -  8 

C"  2pw 

1  -  2  +  3  -  4  + 

5  -  6  +  7  -  8 

A2u 

Fe  4ps 

N  2pir 

1  +  2  +  3  +  4 

C  2p* 

1  +  2  +  3  +  4 

C*  2pw 

1  +  2  +  3  +  4  + 

5+  6+  7  +  8 

C"  2p* 

1  +  2  +  3  +  4  + 

5  +  6  +  7  +  8 

Bl* 

C  2p» 

1  -  2  +  3  -  4 

C‘  2pir 

1  +  2  -  3-  4  + 

5+  6  -  7  -  8 

C”  2pw 

1  +  2  -  3  -  4  + 

5  +  6  -  7  -  8 

B2u 

N  2pw 

1  -  2  +  3  -  4 

C’  2pw 

1  -  2  -  3  +  4  + 

i  -  6  -  7  -  8 

C"  2pir 

1  -  2  -  3+  4  + 

5  -  6  -  7  -  8 

E« 

F«  «„•  t»y.) 

N  2pv 

2-4,  (1-3) 

C  2p* 

1  +  2  -  3  -  4, 

(1  -  2  -  3  +  4) 

C'  2pw 

2  +  3  -  6  -  7. 

(1  -  4  -  5  +  8) 

C’  2pw 

1  +  4  -  5  -  8, 

(2  -  3  -  6  +  7) 

C"  2px 

2  +  3  -  6  -  7, 

(1  -  4  -  5  +  8) 

C"  2pw 

1  +  4  -  5  -  8. 

(2  -  3  -  6  +  7) 

Here  C'  2pir  1-2+3-4+5-6+7  -8,  for  example,  means 


*HAlu)  s  N  U2P»)c.j  *  (zPff)c.2  +  (2P»)c<3  *  (2Pw)C'4  +  (2Pw)c*5 
’  (2Pff)c,6+  ^2pTr^C’7  "  (2pir)c’8  )  ’ 


with  a  suitable  normalisation  constant  N 
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(ll)  Orbitals  of  the  <r- character 


Alg 

Fe 

Fe 

4s 

N 

V 

1  +  2  +  3  +  4 

Blg 

Fe 

3dx* _y» 

N 

V 

1  -  2  +  3  -  4 

B2g 

Fe 

**xy 

Eu 

Fe 

4px,  (4py) 

N 

2-4.  (i-3) 

APPENDIX  H 

Matrix  elements  of  the  potentials  due  to  a 
point  charge,  a  dipole,  and  a  quadrupole 

We  denote  the  radial  part  of  an  atomic  orbital  by  “hi  <r>  •  and 
introduce  the  following  integrals: 

(k)u'  “  it  \  CR«*  R*r  1  r/*  r* J  r  » 

*  |  V,  R,.fCR/r)k"r>ar  ) 

"'R 


for  the  point  charge  case. 


*  k  f  R„l  R..t-  <  R/r  )k"  r1  d  r  1 

'R  * 


(Al) 


(A2) 


for  tha  dipole  case. 


By  uaing  them,  the  matrix  elements  can  be  expressed  as  follows: 

(  A  I  VI  A  )  *  C  0  )AAt 

(  X  IVI  x)  -  (  0)pp  +  1 

(  X  1  VI  x  )  **  ( y  \  v  l  y  )  *  Cojpp  -  yl2Jpp, 

(  l1!  VI  i*)  «  (  0  Jdd  +  —  l  2  J  dd  «-  t  4  j  dd, 

(x*-/l  VU*-/1)  *  Ixy  |  V|Ay  )  -  t  oj  U  -  t  4J  Jj, 

(  y  x  l  V  l  y  a)  »(ix|\/lix)»(o)dj  i-  -i-  i  2.  J  d  j  -  A-  1 4- J  aJ, 

(  a  I  VI  l2jAi  , 

l  .4  IVI  X>  -  JJ*  C  I  ]ip  , 

(  *  I  VI  x»j  -  ^  l  Upd  *  ^  t  33pJ, 

(xlVlixl-(ytVliy)j  l  l  )pd  -  ^  jlL=-  (  3  3  pd . 

In  the  case  of  a  quadrupole,  we  use 

<k)„.  “  -(zkT  •>  R.i  R,.f  1  r/R)kr*Jr, 

tad 

<k,tf  -  |i|  j\,Rvl.< 'VR^i-Mr 

+  R.i  Ru'i’  (  R/r)  r*dr  |  . 


(A3) 


(A4) 


The  matrix  elements  are  as  follows: 


(  A  I  VI  A)  -  (Ojas, 

(  1  I  VI  Z)  -  C  0)fr  4  1.  t23rf  -  -§•  C2)rr  , 

(  x  I  V  I  X  )  -  (  0)  ff  “‘5’C2^fP',’5f  <4?fr  / 

(  ylVI y  )  -  (o3rr  -  ^(j)rr  -  •J'^rp, 

(  x‘lvl**>  -  (o)Al  +%i*7U  -ff  <4><U, 

(  x*-y*l  VU*-y*;  «  (  o)^  -  A.  (4)44  +  ^-(DdJ  -  A- (4)^4, 

(  xy  I  VI  xy)  -  (  x‘-y‘l  VI  x*-y»)  . 

(  yxlVI^*)"*CoJd4  4  t  z-Jdd  —  ^  (4)J4  -  1 2)dJ  — 

l  *  y  |  V  I  *  x)  »  C  0)44  +  “q*  (  2.J44  -  -£p  (  43dJ  +  (a)d<* 

(  x*-y*  |  V I  l*>  -  -  (*)44  +  -*4pC:(4)d<l, 

(  4  I  VI  *)  -  <*3 a4  -  j^r  ( 2 >aJL , 

(  A  I  V  I  x*-y*)  *  Ji 

(  a  I  V  I  **  >  **  jj*  *  , 

(  i  I  V  I  i*3  =  J-y-  C  1  Jr4  +•  C  3 4  -  yj|r  *  3V*, 

(  *  1  Vl^-y)-  o;^ 

( x l Vlix )»  jf-  ClJH  "  t$F  t,,rJ  *$t  <3’H/ 

(y  IVlyiJ-j^-tUyJ  -  C37fJ  -  (  3)fJ. 


-a?- 


The  integrals  [k]^,  and  (k )it,  can  be  expressed  in 
terms  of  elementary  functions  when  the  atomic  orbitals  involved  are  Slater 
functions  and  their  principal  quantum  numbers  are  integers.  When  the 
principal  quantum  number  is  not  an  integer  (e.g.  for  4s  and  4p),  the  inte¬ 
grals  can  be  expressed  in  terms  of  incomplete  gamma  functions. 
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Figure  2  The  co -ordinate  ayetem  and 
numbering  of  atoms. 
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